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Wc study two-body B-meson decays to a charmonium state (r/c, ric{2S) or he) and a K'^ or 
K*^\S92) meson using a sample of 349 fb~^ of data collected with the BABAR detector at the 
PEP-II asymmetric-energy B Factory at SLAG. We measure B{B° ri^K*") = (5.7 ± 0.6(stat) ± 
0.9(syst)) X IQ-^, B(B° ^ rje{2S)K*") < 3.9 x 10"^, B{B+ -> hcK+) x B{hc r/^j) < 4.8 x 10"^ 
and B{B° hcK*°) x B{hc ??c7) < 2.2 x 10"'' at the 90% C.L., and B(??c(25') ^ KKn) = 
(1.9 ± 0.4(stat) ± l.l(syst))%. We also measure the mass and width of the rjc meson to be m{r]c) = 
(2985.8 ± 1.5(stat) ± 3.1(syst)) MeV/c^ and r{ric) = (36.313:6 (stat) ± 4.4(syst)) MeV. 
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PACS numbers: 13.25.Gv, 13.25.Hw 



In the simplest approximation, B decays to a charmo- 
nium state and a if or K* meson arise from the quark- 
level process b ccs and have been observed to occur 
with large rates However several decay modes are 
still poorly known, particularly in the case of singlet 
states such as 77c and he- A better knowledge of the 
relative abundances of the decay to the various charmo- 
nium states allows a deeper understanding of the under- 
lying strong processes and tests of the predictions of mod- 
els such as non-relativistic QCD [2|. In non-relativistic 
QCD, the B decay rates to all P-wave states of charmo- 
nium, Xcj (J — 0,1,2) and he, do not vanish and are 
foreseen to be comparable in magnitude. 

In this document, we study i?- meson decays to 
{KKtt)K+, {KKtt)K*°, t)c1K+ and ricjK*°, from which 
we measure the branching fractions for the following de- 
cay modes: ^ 7?ci^*°, ^ 7]ci2S)K*", B° 
hcK*°, B+ -> hcK+ 3, and r]c{2S) -> KKn. We also 
measure the mass and width of the ijc meson. The he 
meson has recently been discovered by the CLEO Col- 
laboration as a narrow peak at 3524.4 ± 0.7 MeV/c^ in 
the ?7c7 invariant mass distribution in ip(2S) —> rjc'jTT^ de- 
cays 4J , and this observation was confirmed by the E835 
Collaboration [5i]- The ric{2S) state was discovered by the 
Belle Collaboration in B decays to {K'^K^tt^)K 6], and 
subsequently observed in the processes 77 rjci'iS) — > 
KO/iT^Tr^F and e+e" J/V'r?c(25); its mass is (3637 ± 
4) MeV/c^ and its width is (14 ±7) MeV [l|. No branch- 
ing fraction for any 77c (25*) decay mode is yet listed 
by the Particle Data Group [ij. Using the measured 
value for r(r/c(2S') -> 77) x B{r)c{2S) KKn) 0, 
a measurement of B{r]c{2S) — > KKtt) can be used as 
an input to derive r(77c(25) 77), a quantity cal- 
culable in a theoretically clean way within the con- 
ventional framework of QCD: calculations that assume 
B{t]c{2S) KKtt) = B{ric KKtt) lead to values of 
T{ric{2S) 77) smaller than expectations, pointing to 
a possible anomaly in the rjc{2S) decay Q. The branch- 
ing fraction of B^ — > rjcK*^ is currently known with a 
40% uncertainty, (1.2 ± 0.5) x lO'^ [|, while B decays 
to ric{2S)K* and h^K*^*^ have never been observed. The 
Belle Collaboration studied the decay _B+ hcK^ with 
he Vol and reported B{B+ r]c^K+) < 3.8 x 10"^ at 
the 90% C.L. for an invariant mass of the r]c^ pair in the 
range [3.47,3.57] GeV/c^ This limit is comparable to 
analogous limits for Xc2 but significantly smaller than the 
measured branching fractions for B decays to rjc, J/ip , Xco 
or Xci and a kaon [1]. No other or decay modes 
with he have yet been studied. The mass and width of 
the rjc are important parameters in models of the charmo- 
nium spectrum [11.] : the hyperfine separation (?7c, J/i^) is 
directly related to the spin-spin interaction. The r]c mass 
and width measurements reported in the literature [ij 



are often in poor agreement with one another. The listed 
world average for the mass is (2979. 8± 1.2) MeV/c^, with 
measurements ranging from 2969 to 2984 MeV/c^, and 
for the width it is (26.5 ± 3.5) MeV with values ranging 
from 7 to 48 MeV. 

In this analysis we reconstruct the r/c and rjc{2S) in the 
K^K^TT^^ and K^K~tt^ decay modes, the he in its decay 
to 77c7, the in the mode K^^ 



and the if*° in 



K 



*o 



The K^^K^TT^ and K+K-jt^ final states 



are chosen because they are among the easier rjc decay 
modes to reconstruct and have a rather large branching 
fraction, B{ric ~* KKtt) = (7.0 ± 1.2)% 1]. For the 
r]c{2S), the K'^K^tt^ mode is the only decay observed 
so far. The rjc'y decay of the he is chosen because it is 
expected to comprise about half of the total he decay 
width For decays with rjc and he, we measure ratios 
of branching fractions with respect to B{B^ ijcK^) = 
(9.1 ± 1.3)2< lO""* [H, to cancel the 17% uncertainty on 
B{ric KKtt). Similarly, we measure the ratio B{B'^ 
ric{2S)K*°)/B{B+ r?c(2S')is:+), to cancel the unknown 
branching fraction of ric{2S) KKtt. 

The data used in this analysis were collected with the 
BaBAR detector at the PEP-II e+e~ storage rings, and 
correspond to 349 fb^^ of integrated luminosity collected 
at the T{AS) resonance, comprising 384 million BB pairs. 
The BaBAR detector is described elsewhere ^12]. We 
make use of Monte Carlo (MC) simulations based on 
GEANT4 [13. 

The event selection is optimized by maximizing the 
quantity Ns / V Ns + Nb , where Ns {Nb ) represents the 
number of signal (background) candidates surviving the 
selection. Ns is estimated from samples of simulated 
events of B ricK^*\ric KKn decays for B 
{KKTr)K^*\ and B hcK^*\hc rjcl, rjc KKtt for 
B t]cJK'^*\ Nb is estimated from signal sidebands on 
data, defined by the signal candidates with reconstructed 
e+e~ center-of-mass energy farther than 3 standard de- 
viations from the expectation in e~^e~ collisions at the 
Y{4S) peak. Simulated signal events and data are nor- 
malized to each other using the available measurements 
for B decays to rjc and assuming B{B hcK'^*^) = 
1 X 10^5. 

We select events with BB pairs by requiring at least 
four charged tracks, the ratio of the second to the zeroth 
order Fox- Wolfram moment IJ] to be less than 0.2, and 
the total energy of all the charged and neutral particles 
to be greater than 4.5 GeV. 

Charged pion and kaon candidates are reconstructed 
tracks having at least 12 hits in the drift chamber, a 
transverse momentum with respect to the beam direc- 
tion larger than 100 MeV/c, and a distance of closest ap- 
proach to the beam spot smaller than 1.5 cm in the plane 
transverse to the beam axis and 10 cm along the beam 
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axis. We use particle identification provided by measure- 
ments of the energy loss in the tracking devices and the 
Cherenkov detector. A K*'^ candidate is formed from a 
pair of oppositely charged kaon and pion candidates orig- 
inating from a common vertex and having an invariant 
mass within 60 MeV/c^ of the nominal if *° mass 

Photon candidates are energy deposits in the elec- 
tromagnetic calorimeter that are not associated with 
charged tracks, having energy greater than 100 MeV 
and a shower shape consistent with that of a photon. A 
7^0 candidate is formed from a pair of photon candi- 

dates with invariant mass in the range [115,150] MeV/c^ 
and energy greater than 400 MeV. These candidates are 
constrained to the nominal 7r° mass |1| . 

A K'^ — > TT+TT^ candidate is formed from a pair of op- 
positely charged tracks originating from a common vertex 
and having an invariant mass within 20 MeV/ of the K'^ 
mass. Its measured decay-length significance is required 
to exceed three standard deviations. The candidate is 
constrained to the nominal if" mass \X!\ ■ 

The B+'° {KKtt)K+'*" candidates are formed by 
pairing a if*" or candidate, referred to as the pri- 
mary kaon, and a if^if^Tr^ or if+if ~7r" combination 
with invariant mass above 2.75 GeV/c^ to include the 
whole charmonium region. The i?"*"'" rjcjK'^'*'^ candi- 
dates are formed by combining a if*" or if"*" candidate, a 
photon with energy exceeding 250 MeV, and a if "if ^tt^ 
or if """if ~7r" combination with invariant mass consistent 
with the r]c mass. We perform a vertex fit to the B can- 
didates and require the probability to exceed 0.002. 
We define two kinematic variables: the beam-energy sub- 
stituted mass, mEs = \/ ^bcam ~ P% = Eb — 
-Ebeam, whcrc pB {Eb) IS the reconstructed B momen- 
tum (energy) and ii^beam is the beam energy, in the e+e" 
center-of-mass (cm.) frame. B candidates are retained 
if they have mES greater than 5.2 GeV/c^ and Ai? within 
[-24,30], [-40,30], [-34,30], and [-40,30] MeV for the 
K°K^Ti^K*°-+, ii:+A'-7r"if*"'+, K°^K^TT^-fK*°-'+, and 
if +if ~7r"7if *"'+ combinations, respectively. B mesons 
produced in the process T{AS) — > BB follow a sin^ 9b 
distribution, where is the polar angle of the B can- 
didate momentum vector in the e^e'^ cm. frame: we 
require [cos 0^1 < 0.9. 

To suppress background, if+7r^, K^K^, K'^K^, 
if"7r^ and if+7r~7r+ combinations with invariant masses 
within 30 MeV/ of the £)", Dg and meson masses [l| 
are excluded when forming B candidates. We also re- 
move pairs containing a primary kaon where the 
invariant mass of the pair is within 30 MeV/ of the (f) 
meson mass [H. 

In events where more than one B candidate survives 
the selection, the one with the smallest |Ai?| is retained. 
In cases of multiple B candidates composed from the 
same final state particles, and thus having the same value 
of [Ai^l, we retain the one for which the primary kaon 
has the largest momentum in the e^e~ cm. frame. 



The samples surviving the selection include a sig- 
nal component, a combinatorial background component 
given by random combinations of tracks and neutral clus- 
ters both from BB and continuum events e^e~ — s- qq 
{q = u,d,s,c), and a component due to B decays with 
a similar final state to the signal. As opposed to the 
combinatorial background, such "peaking backgrounds" 
exhibit the same distribution as signal events in toes 
and AE, but their ifif7r(7) invariant-mass distribution 
(mx) is different. The signal content in data is there- 
fore obtained by means of a maximum likelihood fit to 
nix for all candidates having toes in the signal region 
[5.274, 5.284] GeV/c^, after subtracting the combinatorial 
background. The mx distribution for the combinatorial 
background events is obtained by extrapolating into the 
Toes signal region the mx distribution measured in the 
Toes sideband, defined by 5.20 < toes < 5.26 GeV/c^. 
The correlation between mx and toes is found to be 
negligible in the relevant regions. A binned fit is then 
performed on the mES-sideband-subtracted mx distri- 
bution. 

To estimate the background we perform an unbinned 
maximum likelihood fit to the toes distribution as fol- 
lows. The B component, accounting for the sum of signal 
and peaking background, is modelled by a Gaussian func- 
tion whose width is taken from the simulation and whose 
mean is fixed to the i3-meson mass [if. The toes dis- 
tribution of the combinatorial background is represented 
by an ARGUS function \15^ . The total number of events 
and the exponent of the ARGUS function are left free 
in the fit. The spectrum for candidates in the toes side- 
band is normalized to the toes signal window by using 
the integrals of the ARGUS component in the two regions 
(Fig.©. 

The mx distribution for B+ {KKtt)K+ and B" 
(if if 7r)if *" is shown in Fig. [2l after subtraction of the 
Toes sideband background. The two samples are simulta- 
neously fitted to the sum of an rjc, an ric{2S), a J/ij} , a Xci 
and a 'iJj{2S), and a background component accounted for 
by first-order polynomials. The r]c and ric{2S) peaks are 
modelled by a non-relativistic Breit-Wigner convolved 
with a Gaussian function, the others by Gaussians. The 
masses of Xcii ??c(25') and ip{2S) and the width of the 
r]c{2S) are fixed to the world average values [H. To re- 
duce systematic uncertainties on the r/c mass measure- 
ment from potential distortion effects in data shifting 
the peak positions, in the fit we float the mass of the 
J/ip and fit for the mass difference between J/ip and rjc- 
We also fioat the width of the rjc, the mass resolutions, 
modelled by the widths of the Gaussian functions, sepa- 
rately for the if"if ='=7r^ and if +if ~7r" modes; the coef- 
ficients of the background polynomial functions and the 
number of signal and background events. The fit ex- 
tends over the mx range [2.75,3.95] GeV/c^. No com- 
ponent is included for other charmonium states such as 
XcOj he and Xc2, since they have not been observed to 
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FIG. 1: The rrtES distributions for (a) B+ {KKn)K+, (b) 
B° {KKit)K*°, (c) B+ r,c'yK+ and (d) 5° rja-yK'^ can- 
didates; points with error bars are data, the soUd hne represents 
the result of the fit described in the text, and the dotted line rep- 
resents the ARGUS background parameterization. No appreciable 
B component, either signal or peaking background, is observed for 
the 5+ ric'yK+ and rjc^K*^ cases. 



decay to KKtt and/or in B decays. Table H] summarizes 
the numbers of events found by the fit, separately for the 
B+ {KKtt)K+ and {KKtt)K*° samples. The 

of the fit divided by the number of degrees of freedom 
{Noof) is 1.2. The mass resolution is determined by the 
fit to be 9 ± 1 MeV/c2 and 20 ± 9 MeV/c^ for K^K^tt^^ 
and K'^K~7t'^, respectively. The mass of the J/ip is 
found to be 3096.4 ± 1.0 MeV/c^, the mass difference be- 
tween J/ip and 7yclll-l±l-5 MeV/c^, and the t]c width 
Se.StH MeV. Using to(J/V') = 3096.916±0.011 MeV/c^ 
from if, we derive 771(77^) = 2985.8 ± 1.5 MeV/c^. 



TABLE I: Numbers of r?c, J/ip , Xa, VcCiS) and 1^(25') events 
obtained from the fit described in the text with statistical 
errors only. 





B+ {KK-!v)K+ 


B" ^ (KKn)K'-'" 




732 ± 27 


189 ± 18 




154 ± 15 


56 ±9 




59± 10 


13 ±7 




59± 12 


13 ±9 




15 ±8 


0±4 



In the case of rjcl^^ ^^^d ijcjK*^ 

(Fig. [3]), the mEs-sideband-subtracted mx distribution is 
fitted to the sum of an he signal modelled by a Gaussian, 
and a background represented by a first-order polyno- 
mial. The mass of the he is fixed to the CLEO measure- 
ment, 3524 MeV/c^ The Gaussian resolution is fixed 




3.2 3.4 3.6 3.8 

mx(GeV/c2) 




3.2 3.4 3.6 3.8 
mx(GeV/c2) 

FIG. 2: Fit result (solid line) superimposed on the rriES-sideband- 
subtracted mx distribution (points with error bars) for (a) B'^ — ► 
{KKtt)K+ and (b) B° {KKn)K*°. 



to the value determined from MC events, 16 MeV/c^ [la |. 
In the fit, the numbers of signal and background events 
are left free. The fit is performed over the mx range 
[3.3,3.7] GeV/c^. It yields 11 ± 6 and 21 ± 8 ftc candi- 
dates with a x^/NdoF of 41/39 and 42/39 for the B+ 
and B'^ yields, respectively. 

The stability of the fit results is verified for vari- 
ous configurations of the fitting conditions. For B 
{KKtt)K^*\ we perform the fits with and without com- 
ponents for XcOj Xc2, he and ip{2S) in various combina- 
tions. The values for the signal yields and the floated 
parameters returned by these fits are consistent with the 
nominal configuration. We validate the fit procedure us- 
ing a MC technique: we simulate a number of exper- 
iments by randomly generating samples of events dis- 
tributed in mx according to the models used in the fit. 
The number of events generated is equal to the number 
of events in the corresponding real data sample. The 
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FIG. 3: Fit result (solid line) superimposed on the mES-sideband- 
subtracted mx distribution (points with error bars) for (a) 
Tjc'yK^ and (b) — > rjc'yK*'^. No significant he signal is evident. 
The dashed line is the result of the fit with no signal component. 



parameters of the distributions are set to their fixed or 
fitted values. The fit is repeated under the same condi- 
tions as used on real data. The numbers of signal and 
background events are distributed as expected. The ro- 
bustness of the fit is tested on simulated events by varying 
the numbers of signal and background events input, in- 
cluding the null result. The numbers of events returned 
by the fit are consistent with the inputs for all cases. As 
additional cross-checks, we verify that the observed num- 
bers of J/ip , Xci and ip{2S) candidates in the data agree 
with the expectations. 

We evaluate systematic uncertainties on the numbers 
of signal candidates and the mass and width determina- 
tion by individually varying the parameters that are fixed 
in the fits by ±1 standard deviation from their nominal 
values. We also estimate the systematic uncertainties 
that arise from a different choice of binning, fit range, and 



background parameterization. For {KKTr)K~^ 
and -8° {KKtt)K*'^ , where the mass resolutions are 
floated, we estimate an additional systematic uncertainty 
by taking the variations with respect to a fit performed by 
fixing the mass resolutions to the values determined from 
the simulation, 8 MeV/c^ and 19 MeV/c^ for JC^X^tt^f 
and K'^K~Tr^, respectively. The large natural widths 
of the r]c and r]c{'2S) introduce the possibility of inter- 
ference effects with non-resonant B decays to the same 
final state particles. This can modify the mx distribu- 
tion with respect to the one used in the fit. The fit is 
repeated including an interference term between the rjc 
and the background in the fitting functions. The am- 
plitude and phase of the interference term are left free 
in the fit. The variation of the r]c yield with respect to 
the nominal fit is taken as an estimate of the systematic 
error due to neglecting interference effects. A similar ap- 
proach is undertaken for 77c (25'). Summing in quadrature 
all the contributions, the total systematic uncertainty on 
the signal yield determination is 6%, 3%, 25%, 18%, 25% 
and 23% for B+ 'ncK+ , B° VcK*°, B+ ^ hcK+ , 
B° hcK*°, B+ r]ci2S)K+ and 5° r^^{2S)K*° , 
respectively, and the total systematic uncertainties on 
the mass and width are 3.1 MeV/c^ and 4.4 MeV, 
respectively. 

The selection efficiency for 5+ VcK'^ is 6%. The 
ratios of the selection efficiencies with respect to 5+ 
r]cK~^, estimated by using simulated events, are, includ- 
ing systematic uncertainties, 0.64 ± 0.01, 0.51 ± 0.01, 
0.29 ± 0.02, 0.84 ± 0.01 and 0.54 ± 0.01 for B° rjcK*°, 
B+ -> KK+, B° hcK*°, B+ ric{2S)K+ and 
Tlc{2S)K*''\ respectively. Most uncertainties on 
the efficiencies cancel out in the ratios because of the sim- 
ilar final states. The remaining uncertainties are mainly 
due to differences between real data and simulation in the 
photon reconstruction as estimated from photon control 
samples from data (1.8%), and the unknown polarization 
for B° hcK*° estimated as in [H] (6%). 

As a check, using the signal efficiency computed from 
MC events, the signal yield observed in data, and the 
number of BB pairs in the data sample, we derive 
B{B+ T]cK+) X B{t]c KKtt) = (8.0 ± 0.4(stat)) x 
10^^. This is in agreement with the world average value 
of (6.4 ± 1.4) X lO-'^ 

We calculate the ratios of the branching fractions with 
respect to B{B^ rjcK^) using the ratios of signal 
yields and efhciencies with respect to i?+ rjcK'^, 
Rr = r(r(4S') ^ B+B-)/V{T{AS) B^B°) = 1.026 ± 
0.032 [1], and B{K*° K+tt-) = 2/3, and summing the 
uncertainties in quadrature. We define Rri^K* = B{B^ — > 
il,K*^)/B{B+ ^ 7^,K+), Rh^K = B{B+ ^ h,K+) x 
B{K -> Vcl)/B{B+ -> r^,K+), Rh^k- - ^(5° ^ 
KK*°) X B{K ^ Vc7)/B{B+ ^ r^,K+)_, i?„,(2S)A' - 
B{B+ ^ Tjc{2S)K+) X B{7]ci2S) ^ KKt:)/{B{B+ ^ 
TjcK+) X Bir]c ^ KKtt)) and R^^i2S)K' = B{B° ^ 
ric{2S)K*°)/B{B+ ric{2S)K+). Table |TT] summarizes 
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the systematic uncertainties on the measurements. 
We obtain R^^k- = 0.62 ± 0.06(stat) ± 0.05(syst), 



R 



u.uyu_Q 



(stat) ± 0.025(syst) and the 90% 
C.L. upper Umits Rh^k < 0.052, Rh^K- < 0.236, and 
Rria(2S}K' < l-O- These are determined by assuming 
that each measurement fonows a Gaussian distribution 
around the central value, with standard deviation given 
by the statistical and systematic uncertainties added in 
quadrature. 

TABLE 11: Summary of the relative contributions to the 
systematic errors on R,,^k' , Rh^K, Rh^K 



R 



Rr]a(2S)K and 



77c{2S)K*- 





O" 


{R)/R (%) 








RrjaK 


Rh^K 


Rh^K' 


Rva{2S)K 


Rria{2S)K* 


Signal yield 


6.6 


26 


19 


26 


34 


Signal efficiency 


1.4 


2.2 


6.7 


1.3 


2.2 


Rr 


3.1 




3.1 




3.1 


Total 


7.2 


26 


20 


26 


34 



Using B{B+ r^cK+) = (9.1 ± 1.3) x lO"'^, we derive 

B{B^ TjcK*°) = (5.7±0.6(stat)±0.4(syst)±0.8(br))xl0' 

where the last error is from the uncertainty on B{B^ 
ijcK'^), and the 90% C.L. upper limits 
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B{B+ 
BiB° 



> hcK+) X B{hc 
hcK*") X B{hc 



rjcj) < 4.8 X 10"^ 
T]cj) < 2.2 X lO""*. 

Vc{2S)K+) 



Using the world average value B{B 
(3.4 ± 1.8) X 10-'';i|, we derive 

B{B" T]c{2S)K*°{890)) < 3.9 x 10" 

at the 90% C.L. Finally, using B{B+ _^ r]cK+) x B{7jc 
KKtt) = (6.88 ± 0.77 " " 

B{ric{2S) ^ KKtt) 



toil) X 10"^ tl8|. we derive 
(1.9±0.4(stat)±0.5(syst)±1.0(br))%, 



where the last error accounts for the uncertainties on the 
branching fractions used in the calculation. 

In summary, we obtain a measurement of B{B'^ — > 
rjcK*^) in agreement with, and CTeatly improving upon, 
the previous world average value [l| . We obtain an upper 
limit for B{B^ hcK^) x Bjhc — > ??c7) in agreement 
with the previous Belle result [10|, and set the first upper 
limit on B{B° hcK*°) x B{hc ?7c7): these confirm 
suppression of he production in B decays. We report the 
first upper hmit on B{B° r]c{2S)K*°) and the first 
measurement of B{j]c{2S) — > KKtt). The latter branch- 
ing fraction is smaller than the corresponding branching 
fraction for rjc, and can be used to derive T{ric{'2S) 77). 
We measure m{ric) = 2985.8 ± 1.5 ± 3.1 MeV/c^ and 
r{r]c) =36.3^16^4.4 MeV. These are in agreement with 
previous BaBAR measurements from 77 collisions [l9j and 
slightly higher than the world average values [l[. 
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